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1. Introduction 

The determination of the equation of state (EOS) of strongly interacting matter on the lattice 
requires actions with small discretization effects, a realistic light quark spectrum and good control 
over the systematic errors in order for the result to be relevant in helping to understand experimental 
findings. For our thermodynamics studies we use the asqtad quark action §\§ for 2+1 flavors, 
combined with a one-loop Symanzik improved gauge action [|fl. Both actions are highly improved 
and have discretization errors of 0(a s a 2 , a 4 ) and 0(a 2 a 2 , a 4 ), respectively. We do our simulations 
along trajectories of constant physics using the dynamical R algorithm [Q] and thus our calculations 
are subject to finite step-size errors. Here we describe the method we use to estimate and correct 
for this systematic error. 

To approximate the experimental conditions as closely as possible we have started a project 
which extends our EOS calculation to a small non-zero chemical potential. For this purpose we use 
the Bielefeld-Swansea Taylor expansion method [Q] and preliminary results for the EOS at a set of 
chemical potentials are presented. 

2. Simulation overview 

We study the quark-gluon system at a range of temperatures along trajectories of constant 
physics. Such trajectories are defined by keeping the ratios m K /m p and m^Jm^ fixed while chang- 
ing a at a constant N t . We work with two approximate trajectories; for both of them the heavy 
(strange) quark mass m s is fixed to the physical value within about 20%. The light quark masses 
for the two trajectories are m lu / m 0.2m s (rn n jm p m 0.4) and m u d m 0.lm s (m^/mp m 0.3), respec- 
tively. For both trajectories we have performed calculations at N t = 6. In addition, a study at N t = 4 
was done for the w O.lm^ trajectory in order to compare the effect of the larger discretization 
errors on the EOS. The constant physics trajectories are parameterized with RG-inspired formulae 
using the hadron spectrum data at certain anchor points (see [^J for explicit formulae). 

The lattice gauge configurations are generated with step sizes in the R algorithm chosen to be 
the smaller of 0.02 and 2m u d/3, and occasionally smaller still. Most of the gauge configurations 
used for the analysis were generated before the RHMC algorithm [§] became known. Since Lattice 
2005 [§], we have doubled the statistics and added a new high-temperature run at j8 = 7.08 (a m 
0.086 fm) and a zero-temperature run at j3 = 6.275 (a rs 0.232 fm) for the N t = 4 case. For both 
trajectories new zero- and high-temperature runs at j8 = 6.85 (a « 0.1 10 fm) were included as well. 
For the list of the rest of the run parameters see [||]. The lattice spacing a has been determined at 
zero temperature using the heavy quark potential ("n") with the overall scale set by the T 1S-2S 
mass splitting [[7p. We fit all the available data to an appropriate one-loop RG-inspired formula [J|] 
which gives the lattice spacing as a function of the quark masses and the gauge coupling and allows 
us to determine the temperature along the constant physics trajectories. 

3. The EOS analytic form 

We determine the EOS using the integral method [Q], where the pressure is calculated as an 
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integral of the interaction measure and the energy density is a linear combination of both: 

la 4 = -6*A (P) - 12#L A {R) _ 16 ^ A (C ) 
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(3.1) 



pa 4 = - r° * I{a'){a') A dmd ', (3.2) 

Jlnao 

ea 4 = Ia 4 + 3pa 4 . (3.3) 
For physics quantities and parameter definitions in the above see 

4. Systematic errors and EOS results 

Our calculation is affected by the following systematic errors: finite volume effects, choice 



of the lower integration limit in Eq. (3.2) and the finite step-size error in the R algorithm. The 
first two errors are relatively straightforward to estimate. By conducting an EOS calculation along 
the m M j ?a 0Am s , N t =4 trajectory on a small (8 3 x 4) volume and comparing it with the result of 



the calculation on a larger (12 x 4) volume, we have determined that the finite volume effects are 



negligible. To estimate the error introduced by postulating the lower integration limit in Eq. ( p.2| ) to 
be at the coarsest available lattice scale in our simulations on a given trajectory, we have calculated 
the pressure of an ideal pion gas at the corresponding lowest available temperature points. The 
value of the ideal pion gas pressure at these temperatures is about as large as the statistical errors 
we have. Thus, we have ignored this systematic error as well. 

The finite step-size error determination is a much more involved procedure. We have per- 
formed some additional simulations at different step sizes using the R algorithm and, recently, 
some using the exact RHMC algorithm. We measured most of the gluonic and the fermionic 
observables in Eq. ( |3~l| ) at different step sizes. Figure [T] shows the step-size dependence of the 
plaquette and the light chiral condensate for one case. The step-size error in the plaquette variable 
has a potentially significant effect on the EOS, whereas the step-size error in the chiral condensate 
(and the rest of the fermionic observables) is negligible. We have found that the gluonic observ- 
ables (plaquette, rectangle and parallelogram) have similar step-size dependence slopes. For this 
reason we have concentrated on studying the finite step-size corrections mainly for the plaquettes 
and applied the same corrections to the rest of the gluonic observables. Figure || shows the pla- 
quette slopes determined from a set of runs with different step sizes for both trajectories. We see 
that for the m U( i k, 0.2m s trajectory the slopes are small which means that the finite step-size errors 
can be ignored in this case. For the m U( j k, 0.lm s trajectory we have more data for the plaquette 
slopes from both N t = 6 and 4 simulations. We fit our data to a quadratic form and use it to cor- 
rect the EOS. Still, even in this case the finite step-size errors are quite small. The finite step-size 
corrections which we showed at the Lattice 2006 conference were overestimated due to the limited 
amount of data on step-size dependence we had at that time. Figure |3| and || show our results for 
the interaction measure, pressure and energy density after the small finite step-size corrections in 
the m u d ?s 0.\m s trajectory have been included. 
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Figure 1: Plaquette (left panel) and chiral condensate (right panel) vs the squared step size e 2 for the 
improved action for the ensemble at j3 = 6.467, am m \ = 0.01676 and am s = 0.0821. The squared step size 
used for production of this ensemble is 0.0001. 
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Figure 2: Plaquette slopes for the m uc / w Q.lm s trajectory (left) and the m Ilc j « Q.2m s trajectory (right) vs. 
the chiral condensate. 



5. The EOS with non-zero chemical potential 

To include a non-zero chemical potential in the EOS calculation we use the Bielefeld 
Taylor expansion method [BJ|. According to this method the pressure can be expanded as 

p__ x j^l_ v ( ^(h\ n (h\ m 

T 4 ~VT^- n ^ =0 Cnm[I) {Tj \T J ' 

where Z is the partition function, and jlu, are the chemical potentials for the light and heavy quarks, 
respectively. We note that the problems of rooted staggered fermions at non-zero chemical potential 
[^] are not relevant here since all the expansion coefficients are evaluated in the /I/ /, = theory. 
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Figure 3: Interaction measure vs. temperature for both constant physics trajectories and iVj's. 
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Figure 4: Energy density (left) and pressure (right) vs. temperature for both constant physics trajectories 
and N t 's. 

Due to CP symmetry the non-zero terms in the series have even n + m. The non-zero coefficients 
in the above are 

1 1 d n+m lnZ 



n\ ml Nf d(n,N,) n d(n h N t ) r ' 



(5.2) 



JU/,A=0 



where now the jj^/j, are the chemical potentials in lattice units. Similarly, for the interaction measure 
we have: 

/ _ NfdXnZ 

f4 ~ 



Yj b nm(T) 



Nf din a 

where again only terms with even n + m are non-zero and 



(5.3) 
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(5.4) 



To determine the EOS, we need to calculate derivatives of the asqtad fermion matrix such as 
d n lndetM Ul 



(5.5) 
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These derivatives are estimated on the ensembles of lattices along a constant physics trajectory us- 
ing 200 random sources in the region of the phase transition/crossover and 100 sources outside that 
region. We intend to study the EOS up to 0(jU 6 ) on both constant physics trajectories. Currently we 
have data for only one of them and not enough statistics to resolve the sixth-order terms. Thus the 
preliminary results for the EOS we present here are calculated to 0(/i 4 ) on the m uc \ w 0.lm s , N t = 4 
trajectory only. We also set \ih = here. Figure |5] shows our result for some of the coefficients 



involved in Eq. ( |5.1[ ). An interesting feature is that the coefficients quickly reach the continuum 
Stefan-Boltzmann limit above T c . The pressure and energy density corrections due to the non-zero 
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Figure 5: Some Taylor coefficients in the pressure expansion for the m ui j w 0. lm s , N t — 4 trajectory, 
chemical potential are shown in Figure H. 



6. Conclusions 

We have calculated the EOS along two constant physics trajectories with m U( j m 0.\m s and 
0.2m. s using the R algorithm. To estimate the finite step-size errors we have performed a number 
of additional R algorithm simulations at different step-sizes and a few RHMC simulations along 
the constant physics trajectories. The analysis of this additional data shows that the finite step-size 
errors are negligible on the m u d ~ 0.2m. s trajectory and small (in most cases less a than few percent) 
on the m uc } w 0.1m,, trajectory. 
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We have started a non-zero chemical potential study of the EOS with 2+1 flavors using the 
Bielefeld-Swansea Taylor expansion method and presented preliminary EOS results for a number 
of different \i\jT values. 
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